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METHOD AND APPARATUS FOR FAST 
NATURAL LOG(X) CALCULATION 

BACKGROUND OF THE INVENTION 

The invention relates generally to methods and apparatus for 
computing a computationally intensive algorithm, and more specifically to a method 
and apparatus for computing log(x), or equivalently, -log(x), in a manner that is 
particularly useful for computed tomographic image processing and other 
applications. 

In at least one known computed tomography (CT) imaging system 
configuration, an x-ray source projects a fan-shaped beam which is collimated to lie 
within an X-Y plane of a Cartesian coordinate system and generally referred to as the 
"imaging plane. 11 The x-ray beam passes through the object being imaged, such as a 
patient. The beam, after being attenuated by the object, impinges upon an array of 
radiation detectors. The intensity of the attenuated beam radiation received at the 
detector array is dependent upon the attenuation of the x-ray beam by the object. Each 
detector element of the array produces a separate electrical signal that is a 
measurement of the beam attenuation at the detector location. The attenuation 
measurements from all the detectors are acquired separately to produce a transmission 
profile. 

In known third generation CT systems, the x-ray source and the 
detector array are rotated with a gantry within the imaging plane and around the object 
to be imaged so that the angle at which the x-ray beam intersects the object constantly 
changes. A group of x-ray attenuation measurements, i.e., projection data, from the 
detector array at one gantry angle is referred to as a "view." A "scan" of the object 
comprises a set of views made at a different gantry angles, or view angles, during one 
revolution of the x-ray source and detector. In an axial scan, the projection data is 
processed to construct an image that corresponds to a two dimensional slice taken 
through the object. One method for reconstructing an image from a set of projection 



15-CT-5271 




data is referred to in the art as the filtered back projection technique. This process 
converts the attenuation measurements from a scan into intergers called "CT numbers" 
or "Hounsfield units," which are used to control the brightness of a corresponding 
pixel on a cathode ray tube display. 

The negative natural logarithm function -log(x) is an important yet 
computationally intensive algorithm in computed tomographic (CT) image processing. 
In known systems, a 5th order polynomial is used to approximate the function. 
However, this polynomial still consumes more than 20% of the total image processing 
time and generates a relatively large approximation error and error standard deviation. 




Using equation (1), -log(x) can be written as: 

y = - log(x) = - log(m) -ex log(2) (2) 

The following equation uses a finite order polynomial to approximate 
log(w) in a region 1 < m < 2 . Generally speaking, the higher the order of the 
polynomial, the better the approximation will be, but the computational load is in 
proportion to the order of the polynomial. For example, a 5th order polynomial 
presently used is written as: 

log(w) « (a 0 + a } m + a 2 m 2 + a 3 m 3 + a 4 m 4 + a 5 m 5 ) (3a) 
or as: 



y = - log(;c) « ~(a 0 + m(a ] + m(a 2 + m(a 3 + m{a A + a $ m)))) + e x a 6 ) 



(3b) 
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In equation (3b), a 0 — a 6 are precalculated constants. To compute 

-log(jc), six additions and six multiplications are required, plus mantissa and exponent 
extractions. 

To process images more efficiently and accurately, it would be 
desirable to provide methods and apparatus to reduce the complexity of the 
approximation used to calculate -log(:c) while achieving numerical accuracy consistent 
with IEEE (Institute of Electrical and Electronic Engineers) floating-point precision. 

BRIEF SUMMARY OF THE INVENTION 




""S^ere is therefore provided, in one embodiment of the present 
invention, a method forcomputing a natural logarithm function that includes steps of: 
partitioning a mantissa re^km between 1 and 2 into N equally spaced sub-regions; 
precomputing centerpoints a, of^ch of the N equally spaced sub-regions, where 
i = 0 5 . . N -1 ; selecting N sufficiently n&ge so that, for each sub-region, a first degree 
polynomial in m computes log(/w) to within a preselected degree of accuracy for any m 
within the sub-region, where m is a mantissa of a binary floating point representation 
of a number; and computing a value of log(x) for a binary Abating point representation 
of a particular number x stored in a memory of a computing device utilizing the first 
degree polynomial in m. 

It will be seen that this embodiment and others described herein reduce 
the complexity of approximations used to calculate natural logarithms while achieving 
numerical accuracy consistent with EEEE floating point precision. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a pictorial view of a CT imaging system. 



Figure 2 is a block schematic diagram of the system illustrated in 

Figure 1 . 




DETAILED DESCRIPTION OF THE INVENTION 

Referring to Figures 1 and 2, a computed tomography (CT) imaging 
system 10 is shown as including a gantry 12 representative of a "third generation" CT 
scanner. Gantry 12 has an x-ray source 14 that projects a beam of x-rays 16 toward a 
detector array 18 on the opposite side of gantry 12. Detector array 18 is formed by 
detector elements 20 which together sense the projected x-rays that pass through an 
object 22, for example a medical patient. Detector array 18 may be fabricated in a 
single slice or multi-slice configuration. Each detector element 20 produces an 
electrical signal that represents the intensity of an impinging x-ray beam and hence the 
attenuation of the beam as it passes through patient 22. During a scan to acquire x-ray 
projection data, gantry 12 and the components mounted thereon rotate about a center 
of rotation 24. 

Rotation of gantry 12 and the operation of x-ray source 14 are 
governed by a control mechanism 26 of CT system 10. Control mechanism 26 
includes an x-ray controller 28 that provides power and timing signals to x-ray source 
14 and a gantry motor controller 30 that controls the rotational speed and position of 
gantry 12. A data acquisition system (DAS) 32 in control mechanism 26 samples 
analog data from detector elements 20 and converts the data to digital signals for 
subsequent processing. An image reconstructor 34 receives sampled and digitized x- 
ray data from DAS 32 and performs high speed image reconstruction. The 
reconstructed image is applied as an input to a computer 36 which stores the image in 
a mass storage device 38. 

Computer 36 also receives commands and scanning parameters from 
an operator via console 40 that has a keyboard. An associated cathode ray tube 
display 42 allows the operator to observe the reconstructed image and other data from 
computer 36. The operator supplied commands and parameters are used by computer 
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36 to provide control signals and information to DAS 32, x-ray controller 28 and 
gantry motor controller 30. In addition, computer 36 operates a table motor controller 
44 which controls a motorized table 46 to position patient 22 in gantry 12. 
Particularly, table 46 moves portions of patient 22 through gantry opening 48. 

A negative natural logarithm function -\og(m) is used by image 
reconstructor 34 to generate images. In one embodiment of the invention, the function 
log(m) of equation (3a) above is written as: 



log(/w) = log(a) + 
log(m) (n) 



+ 



TV! 



1! 2! 



(4a) 



or as: 



10 



i /x , /x (m-a) (m - a) 2 
log(/w) « log(a) + ^ L - ± ^— + 



a 



2a 2 



(4b) 



written as: 



where a is a known reference point. The error of the above function is 




error < 



log(m) 



(«+o| 



N\ 



(m-a) 



n+l 



(5) 



because (m-a)<l, there are two ways to minimize the error. One way 
15 is fo increase the orderbf s the approximation, and the other is to minimize the distance 
from m to a. Because mantiSsa^m is between 1 and 2, in one embodiment of the 
present invention, the region betweenl^aQd 2 is partitioned into N equally spaced sub- 
regions. Centers of each of the sub-regions are^recomputed and used as reference 
points in equations (4a) and (4b). By partitioning infcKa^sufficiently large number of 
20 sub-regions, a low order polynomial function produces sufficient accuracy for CT 
imaging purposes. In particular, by selecting a sufficiently large^number of sub- 
regions, for any m within any particular sub-region, log(w) is comput^dsby a first- 
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de^ee^olyHomial to within a preselected degree of accuracy within that sub-region. 
For example, computer 36 uses the first de^ee-pGlyoomial in m to compute values of 
log(x) for binary floating point representations of particularlltrmhgrs x stored in its 
memory. ^ ^^^^ 

A log(/w) approximation that is based on a first order polynomial with a 
set of precalculated reference points is written as follows: 

log(m)^log(a t )+ ' (/W ~ g|) ; i=0,...,tf-l; l<a { <2 (6) 

<*i 

where a, is a closest reference point to a given mantissa m. 

^Rather than compute a sub-region index using / = round ((m -l)xiV), 
which would require six operations, one embodiment of the present invention reduces 
computation load as follows. A partitioning algorithm divides the mantissa of a 
binary floating point rmmber in memory into two sub-regions. The sub-regions have 
index i and Ac, where Ax ft^a distance from mantissa m to reference point a,-. Indices / 
and Ax are directly extracted from an IEEE floating-point number stored in a computer 
system, thereby reducing computation time and improving accuracy. In one 
embodiment, mantissa partitioning > oGcurs as illustrated in Figure 3, in which index / 
ranges from 0 to 127 and each region repi^sents information extracted from the datum 
shown in Figure 3. More particularly, in^ single precision IEEE floating point 
number, b$i represents a sign bit, b3o the mosk significant bit of exponent e, 6^5 the 
least significant bit of exponent e, b22 the most significant bit of mantissa m 7 and bo 
the least significant bit of mantissa m. (If it is desired to use a different designation 
for the numbering of bits 6, those skilled in the art can make the appropriate changes 
required in the description for notational consistency.) vn this single precision 
embodiment, exponent e is extracted directly from bits b$ 0 to ^'region i is extracted 
directly from bits 622 to bj6\ and Ac (a distance from mantissa m to^eference point a,-) 
is extracted directly from bits bis to bo. 
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UsimT"ThE~"e?cte^ in Figu re 3, a maximum error of 

equation (6) in each sub-region is estimated by an expression written.as: 



error < - 



-^-xf— 1 ; \ = 0,... y N-l; l<a,. <2 
2a, 2 {2NJ \ 



(7a) 



From equation (7a), it is seen that the error of the first order 
approximation is always positive, so that the error is biased. To minimize the 
maximum error, in one embodiment the mean error of equation (7a) is subtracted from 
equation (6). Thus, the unbiased error is written as: 



\error\ < 



1 



4Na : 



i=0,. ...W-l; l<a, <2 



(7b) 
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Subtracting equation (7b) from equation (6) results in an unbiased first 
order polynomial function for -log(x) written as: 



I 

y = - log(x) « b- t + C/Ax + e x log(2) 



for i = 0,...,Af-l 



(8) 



b, = -log(a,) + 

C; =-l/fl,' 



f 1 ^ 



2N)a, 



(9) 



15 



where a = 1 + ' + , and. Ax is a distance from mantissa m to 
N 



20 



reference point a t . The value Ax is extracted directly from an IEEE floating point 
datum. In one embodiment, log(2) and the b t and c, are pre-calculated and saved in a 
look-up table at initialization time. In one embodiment, the values bj are determined 
from precomputed values of log(a,). For purposes of comparison, computation of 
equation (8) requires only 1/3 as much time as is required to calculate the 5th order 
approximation of equation (3b). 
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Taking into account the relationship of m and Ax, an approximation to 
log(;t) (or equivalently -log(jc)) is computed using a polynomial of first degree in m 
and a precomputed value of log(a,-). 

In one embodiment of the present invention, image reconstructor 34 is 
configured with software or firmware to compute logarithms using one or more of the 
methods of the present invention, when CT imaging system 10 images an object from 
acquired projection data. 

From the preceding description of various embodiments of the present 
invention, it is evident that the complexity of the approximation used to calculate 
-log(x) is reduced, while numerical accuracy consistent with IEEE (Institute of 
Electrical and Electronic Engineers) floating-point precision is maintained. Thus, 
images processed by CT imaging system 10 are processed more efficiently, and 
without loss of detail. Although particular embodiments of the invention have been 
described and illustrated in detail, it is to be clearly understood that the same is 
intended by way of illustration and example only and is not to be taken by way of 
limitation. For example, embodiments of the improved computation for log(;c) and 
-log(x) can be incorporated into any computational system requiring increased 
efficiency while maintaining computational accuracy. In addition, the present 
invention is suitable for use with floating point numbers having greater or lesser 
precision than those discussed in detail in this description. The modifications 
necessary to accommodate such different precisions will be apparent to those skilled 
in the art, once the invention described herein is thoroughly understood. Accordingly, 
the spirit and scope of the invention are to be limited only by the terms of the 
appended claims and legal equivalents. 



